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On the Asymptotic Solution of the Non-Homogeneous 

Linear Differential Equation of the n-th 

Order, A Particular Solution, 

By W. Van N. Gaeretson. 



The asymptotic development * for the irregular integrals of a homo- 
geneous linear differential equation has been obtained by both Hornf and 
Love, t Horn has published several papers on the case where the roots of the 
characteristic equation are all distinct, while Love has taken up the case where 
the roots of the characteristic equation are unrestricted as to their order of 
multiplicity, including the case of distinct roots as a special case. 

In this paper we shall consider the non-homogeneous equation where the 
roots of the characteristic equation are distinct, and follow, at the outset, the 
method employed by Dini^ in his researches on linear differential equations. 
In Section I the two theorems stated and proved by him will be generalized so 
as to apply to the non-homogeneous equation and combined in one theorem. 
In Section II we shall determine a particular solution of the given equation. 
To this end we shall make use of the researches of Love $ in the homogeneous 
linear differential equation by employing his solutions. The particular solu- 
tion thus obtained of the non-homogeneous equation will be in the form of 
quadratures. The determination of the asymptotic development of the par- 
ticular integral found in Section II will form the content of Section III. 

Section I. 
Take for consideration the non-homogeneous linear differential equation 

^<»' + %(«;) «/<"-!' + a, (ic) «/<"-2> + . . . .aAco)p=Xix), (1) 

* "Asymptotic Development in Poincare's Sense," of. Acta Mathematica, Vol. VIII (1886), p. 297. 

■If Journal filr Mathematik, Vol. CXXXVIII (1910), pp. 159-191. 

% Annals of Mathematics, Second Series, Vol. XV (1914), pp. 14.5-156. Also American Journal of 
Mathematics, Vol. XXXVI, No. 2 (1914), pp. 151-166. 

lAnnali di Matematioa, Ser. 3, Vol. II (1898), pp. 297-324. /bid., Vol. Ill (1899), pp. 125-183. 
Important contributions have also been made by Poincaire, Kneser, Birkhoff, and others. 
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in which the coefBcients are real or complex functions developable, asymp- 
totically, for large values of x in the form 



(!,•(») ~a?'' 



"[«.o+^^+^+.-..]; i^l,2,....,n; 



k = 0, 1, 2, 



while the first n — i derivatives also possess asymptotic developments. The 
function X{x) will be considered as capable of asymptotic development in the 
form 



X(x)~.-[6.+ ^ + ^ + ....], 



where m, b^, bi, b^, . . . ., are real or complex quantities and bo^O. 

We regard y for the present as a known solution of (1). Let us choose 
n auxiliary functions Si, . . . . , s„ of x, which, with their first n derivatives, are 
continuous for large values of x, and such that for the same values of x the 
determinant 



Q{x). 



Zi Z^ Zi Sj 
^2 ^2 02 • • • • 



.»2 



,(,n-l) 



(2) 



never vanishes. Place 

r[2/"" + «i(a^i)2/"'"" + «2(«i)2/'"~"+ • • • • +«n(a'i)?/]^i(^i)^a;i 

X{x-i)Zi{Xi)dx-i^; i — l,....,n (3) 

where /3^ are sufficiently large, positive, real quantities, or infinity, and x takes 
on large real values.* It is understood that, for the interval of integration, 
all the functions here appearing including y and its first n derivatives are 
continuous. 

An integration by parts gives 

Pi.oy"'-''+Pi.iy''"~''+Pi.2y"'-'' + - ■ ■ ■ +Pi.n-iy 

= f [y{xi)Ziixi)+X{xi)Zi{xj^)]dXi + Ci, (4) 



where 



Pt.O — ^i) Pi.l — ^i^l ^i — ^i% Pi,0> ■•■•> Pi,n-l — ^i'^n-l Pi,n-ii 
-^^ = 0^ — K«-l = «A + fl(^A-l)' + f2(«J«n-2)"+ + Vf>, 

Si={—1Y; i = l, ,n. 



(5) 



*In treating the homogeneous equation, Dini takes /3i = /32=. .. .= /3„, while Horn, in a similar 
discussion, places, as convenience demands, /3i = /32 = . .. .=^/3»=:oo and p„^^ = . .. .=:^„:^^ where 0^v<n 
and /3 is real and positive. The second method is the only one which seems applicable, to the problem in 
hand, in order to insure the convergence of all of the integrals. 
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For the values of x under considei'ation the determinant 

Pl.O Pl.l Pl.n-l 



P{oo) 



never vanishes, because 



Pn.f) Pn.l- 



■Pn.n-1 



P{x)^{-l)'^Q{x). 



The value of y from equations (4) is found to be 

y{x)=g{x)+I, {y{Xi)ki{x,ai)-\-X{Xi)Kiix,Xi))dx, 



(6) 



where 






Si z[ sf-'> Ci 



So. . . . z< 



(n-2) 






/^il^, «i)- ^(.^) and ii:,(x,xi)_ ^^^^ , 

g^(ri;) is the cof actor with respect to the i-th row of the last column of Q{x). 

The formula (6) may be used as a recursion formula. By repeated 
application of the recursion formula the value of y as given in (6) can be 
expressed in the form of an infinite series as follows : 



y{x) = i:^^(x)+x^{x), 



(7) 
where 

n ^x n pXi n /*X\_j^ 

iU^(«)=S ki{x, x^) dx^ • H ki{Xi,x^)dx^ -S Ki<^x~i,^\)ff{^\)<i^\, 

n ^x n ^Xi n p^\—% 

Tx=S ki{x,x^)dxi-ll lii{x^,x^)dXi S A;i(%_2, ^x_i)^a;x_i 

•S Ki{Xy^_^Xy)X{x^)dxj,, ;i = l, 2, 3, 

In the above /Uo(^) =ff{^) and 'To(^) —0. 

The expression for y as given in (7) is still regarded as a known solution 
of (1). The form of the solution being here obtained, we can proceed at once 
to state how an unknown solution could be built up. Choose n auxiliary func- 
tions Si, S2, . . . . , z^ oi X which, with their first n derivatives, are continuous 
for large, positive, real A^alues of x, and such that for the same values of x the 
determinant Q{x) as given in (2) never vanishes. The functions Zi, . . . ., Z„ 
of X can now be formed according to (5); also A{x), g{x), ki{x, x^) and 
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Ki{x, Xi) as given in (6). In building up A{x) it should be noted that 
Ci ,...., c„ are now arbitrary constants and not functions of the betas. To 
simplify the work choose all the c's equal to zero except c^. The function g{x) 

c^A^{x) 



will become 5'r(*) 



Q{x) 



where Ar{x) is the cofactor with respect to c, in 



A{x). Eepresent the corresponding values of /Ux(a;) and y{x) by [Xx,r{^) ^^(^ 
yr{x) respectively. 

We shall now state and prove the following : 

Theorem I. Suppose that a large, positive, real number (i can he found 
such that for the values of (3^ either equal to /? or oo, and for all values of x 
greater than (3 the series 

00 

2/,(a;) = S(^,.,(a;)+T,(a;)) (8) 



x=o 



satisfies the following conditions : 

(a) The series 1iixx,r{'^) and IlTx(a;) converge. 

(b) The series for «/^(«) when multiplied by ]ii{x, x^) may be integrated 

term by term with respect to x-^ from /3 to x, or from a; to oo in case 
f^i equals <x> i = l, . . . .,v where 0<v <7Z. 

(c) The series (8) defines a function y^ix) such that each of the integrals 

y,{x)Zi{x)dx and j Zi{x)X{x)dx (or j yXx')Zi{x)dx and 
2i(a;)Z(a;)(ia/ in case /?i=:oo ; i = \,....,v where ^<v<n) has 

X 

a meaning when x>^. 
Then for such values of x the function yr{x) is an integral of (1). 

Proof: The values of p^^^, i = l, ....,«,- k = l, . . . ., n—1 are given in 
(5). Place 

4>.(«)= r^r(«)-^.(^)^«» s = l, 2, , r—1, r + 1, ,n, 

d>^{x)^ C y,{x)Z,{x)dx-\-c„- 4^,{x)= C z,{x)Xix)dx, s = l, , n, 

Pl,0 Pl,l Pl,n-2 ^l+4'l 



A.(^)' 



A {x) = 



Pn.O Pn,l Pn.n-2 1>n + 4'n 

Pl.O Pl,l Pl,n~l 

Pn.a Pn.X Pn,n~l 



^{-l)'^Q{x). 
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Now by condition (c) the series (8) may be written 

1=1 "^iSi 

By substituting the values of g^{x')'ki{x.,x-^ and JS'j(a;, %) in (9), and taking 
account of equations (3) we find that 

2/,(a;)=A,(a;)/A(a;), 
so it suffices for our proof to show that this function is an integral (1). 

To do this consider the system of n functions >7o > ^yi > • • • • , 57„_i , each 
defined for all values of x sufficiently large by means of the following system 
of n linear equations : 

2',,0>7n-l + 2'»,l>7«-2 + ?',,2'7«-3+ +?'..»-a>7o = ^s + '4's, 5 = 1, 2, , OT. (10) 

It follows that A / \ / A / \ / \ /-1 1 \ 

>7o=A,(a;)/A(a;)=2/,(«). (11) 

By differentiating (10) with respect to x and making use of (5) and (11), we 
find that 

0,(0-Z)+Koei + Kl^2+K2e3+....+PU-20«-l = O, S = l, ....,W, (12) 

where 

^ =>7«-i + «i(«)>7«-2 + «3(«)'7«-3+ +«„-i(aj)>7;+a„(a;)>7o, (13) 

6i=>7«-i— »7«-2; ^i=yin-'i.~n'n~i\ — ; ^n-\—nv—n'^- (14) 

The system (13) consists of n homogeneous equations in n unknowns. 

Upon noting that p^,6 = 2d«i+i— P^.s+iJ^Zq' i' ' "' _2 ^^e discriminant of 

the system reduces at once to ( — \Y~^Q{oi), and hence does not vanish for any 
values of x under consideration, whence 

or, by (11) and (14), 

y!s=y^'\ s=i, 2, — ,n. 

Substituting these values in (13), we find 

y','^^^a^{x)y^"-^> + a,{x)y',-'->+ . . . .+a,{x)y, = X{x), 

which was to be proved. 

Section II. 

In this section we shall endeavor to choose the auxiliary functions Si . . . . s„ 
in such a way that the above theorem may be employed to obtain a particular 
solution of (1). The Complementary function corresponding to (1), i. e., the 
general solution of 

2/"" + %(^)y""-" + «2(«)^'"-"+ • . . • +a,{x)y^Q (15) 
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has already been obtained by Horn and Love. We shall suppose the charac- 
teristic equation of (1), viz., 

w" + ai,oW""-^+a2,o«*"~^+ +««,o = (16) 

to have n distinct roots Wi, ....,%>„. Consider, now, the functions Z^{Xi) 
which occur in kf{x, Xi), ^ = l, . . . ., n. These functions are all of the form of 
the left member of the differential equation 

a„z-a,^,z' + a,_,s"-a„_,z"'+ .... + (-l)"s'«' = 0, (18) 

an equation of the same type as (15). A fundamental system of solutions of 
(18) will be, according to the results obtained by Love,* 

Zi{x)—e-'*'''^x-''-"''Pi_o{^), i — 1, ,n, (19) 

where 

/,.(a;)=M;i^-— ■ +a,-__;i-^ +ai^_k+ij~j^ + +c(.{,-iX, (20) 

and P{_o is of the formf 

Pi,o{^) —(>i.o+ ^~ + "^ + "1 '^ ' ^ — -*-' ''^' 

where p is an arbitrary positive integer and lim 6(,p(«) =0. 

By choosing the auxiliary functions as in (19), viz., the solutions of (18), 
the functions Zj (a;) vanish, and this assures the vanishing of ki{x,Xi). The 
value of «/r(a^) as found in (8) then becomes 

2/r(*)=.9r(a^)+ S C Ki{x,x,)X{Xi)dxi. (21) 

The first term of the right member of (21), viz., g^ix), is a particular solution 
of the homogeneous differential equation (15). To obtain a particular solu- 
tion of (1) set c,. in g^ equal to zero, thus making g^ vanish. Then we have as 
a particular solution of (1), 

n ^x 

y{x)=Yi I K^{x,x^X{x-^dx■^, i = l, ....,n, (22) 

provided the conditions of Theorem I are satisfied. Of these three conditions, 
(a) is satisfied inasmuch as the series reduces to a single term; (b) is satis- 
fied because the functions ^^{x^x-^, % = \, . . . ., n vanish. It will be evident 
that (c) is satisfied when we obtain the asymptotic development for the right 
member of (22). 

* Annals of Mathematics, Second Series, Vol. XV (1914), pp. 145-15G. Also American Joubkal op 
Mathematics, Vol. XXXVI, No. 2 (1914), pp. 1.51-1G6. 

fin this paper the functions P (a)), generally written with subscripts, will have the form of Pi, (,(<») 
here given. 
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By taking the auxiliary functions as above stated, and replacing X{x) by 
its value, the expression for y{<ic) in (22) becomes 

y{x)=ll e^^'-^^x"" "-("-"ffep. ^ (cc) e-^«<*'>a;f -"'^ "P^ Jx^) dx^ . (23 ) 

Section III. 
The Asymptotic Development of y{x). 
Taking?/ (a;) as given in (23), the discussion of its asymptotic develop- 
ment will fall into two parts — Case 1 and Case 2 — according to the behavior 
of /,(^). 

Case 1. 

Suppose /j(»)^0, ^ = l, ....,«. We shall order the f's so that* 
R[h{^)}>R[f2{x)]>B{f^{x)]>....>R[f,{x)}>Q 

>RlUii^)]>R[f.+A^)]>. . ■ .>R[fA^)] 

when X is large, real, and positive. Let |5f, i = l, . . . ., n be so chosen that 
A = /?2= • • • • =f^v = '^ aiid /?^^i = ^^^2= .... =(3^=^, where /3 is to be taken 
sufficiently large. Then «/(«) in (23) becomes 

y{x) = i8,{x)+ h T,{x), (24) 

where 

J, CO 

T^{x) = r%/i<*V<.o-(«-i)*p, ^(a;)e-^''^''a;r-"'-«P^ 2(%)c?a;i, r = v + l, , n. 

Let Sf be the dominant term in 8i, i = l, ...., v. It can be written in the form 
s,{x) =e^«c-V.--<-i)* j .^^ [-/'(^i)]^%, i = l,....,v, 

in which /i(«i)^0, since /,(a3i)^0 for % sufficiently large. An integration 
by parts gives 

s^{x')■= e/*<^'a;"" ->-<«-i'* [/5 ._ i_ oic™-"*- o-*^ 

+ • • • • +A,i.p...*r""°-'-''"'[l+eM(a;i)]]e-^''^'']: 

/y, „m— o,, 0— fc— 1 

a; ' ' 

where 

Pi, 1, • • • • Pt, 1, p„ 1 ) Yi, 1, • • • • ^l, 1, pi, , 

*By R{<x) we shaU mean the real part of a. 
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are easily determined constants, and 

lim 6,_i(a;) =lim Si2{x) =0. 

By repeated application of integration by parts the exponents of the x's 
throughout are continually decreased. After I applications of integration by 
parts, the part still affected by the integral sign, call it Ii{x), is the following: 



^m—ai,(,—Uk+l)—l 
, i, 1^1 



+ •■■• +yu, pw %"'""""""*+"""" [1 + f.-, 2i (%)] )e-''''''^%, i = l,....,v, 
where the y's are easily determined constants and lim ei_2;(^) =0. Suppose I is 

chosen so large that R{m — ai,o~^^ + l))5 — 2. Since e~-^'<*> i = l, . . . ., v is a 
monotonic decreasing function of x, a positive quantity M^ may be found such 

that 

M 

|7,(^.)|<|a;'»-»^|^,,^,;„_„. (27) 

For I greater than unity the last factor here appearing may be made small at 
pleasure with large values of x. By referring to (25), (26) and (27), it is 
readily seen that s^ and consequently Si, i = l, . . . ., v can be represented 
asymptotically in the form 

■U-)~[a^+^ + ^ + ■■■■]■ (28) 

Let ti be the dominant term in T^, i = v + l, . . . ., n. An integration by 
parts gives 

+ ■■■■ +i3,.i,p...a;r""»-'=-'"-[l + f,,i(a;a)] )e-/*<^=']^ 

+ ■■■■ +ni,po,^r""»-*-'""-Hl + e.,.(%)])e-^'^'^'«f*i, i = v + l, ....,n, (29) 

where /^.-^i^o, , A,i,p*,i and yi.i_o, » ni.Pt.i ^^^ easily determined con- 
stants, and, 

lim £{^i{x) =lim Si2{<^) =0. 

After I applications of the integrations by parts, the integral that still remains 
is as follows : 

+ • • • • +yu,Pi, . ^r"""-"'*"-""' [l + e,,2i(«i) ] ) e-^'*''''^«i , 

i — v + 1, ,n. (30) 
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For large values of (3, and for x>^, each term in the integrand of (30) is a 
monotone increasing function of x^ in the interval of integration. A positive 
quantity M^ can be found such that 

M- 

|Z,(a;)|<|^'»-"^| ^„_„;^,, , i = v + l,....,n. (31) 

For I greater than unity the last factor here appearing can be made small at 
pleasure for large values of x. From (29), (30) and (31) it follows that t^, 
and consequently T^, i=y + l, . . . .,n, can be developed asymptotically in the 
form 

T.(.)~[a,.+ ^ + ^ + ....]. (32) 

Case 2. 

Suppose /r(a;)=0 where r is one and only one of the integers 1, . . . ., n. 
Only one / could be identically zero for the roots of the characteristic equa- 
tion (16), are distinct. 

Referring to (23), the part arising from f^ix), in the expression for y{x), 
is as follows : 

u(x) =a>::-'-">[c„,+ ^ +■■■■ + "'■'+;''"' ) 
where p and q are arbitrary positive integers and lim Sp{x) =lim Sg{x)' = 0. It 

x=aa a:=oo 

follows that 

«(a;)~x«-»^+^+^[5o+ § + § + ....], 

provided that m—a^^o^s, when s takes on the values — 1, 0, 1, . . . ., q — 1. 
If, however, m — a^^o^s, then the development of u{x) is 

uix)^x"'-'"'+''+^\Bo+ ^ + 1 + a;"" »-"*=+*= Tco + -~ + "^ + 1 log a;. 

That part of the expansion of (23) arising from i = l, . . . ., r—1, r + 1, ....,% 
takes the same form as (28) or (32) of Case 1. 

In summary we are able to state the following: 

Theoeem II. In the differential equation 

2/'"' + «i(^)^'"~" + «2(^) «/'""''+ • • • • +a„{x)y=X{x), (1) 

45 
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suppose that a^ and X are real or complex functions of x developable asymp- 
totically, when X is large, real, and positive, in the forms 

aAx)~x<'^a,,,+ ^ + ^ + ....~^, i=\,....,n; *=0,1,2,.... 

X(.)~.-[6. +A + _^ ]. ^.*0. 

while the first n — i derivatives of ai{x) also possess asymptotic developments. 
Consider the roots Wi. . . .w^ of the characteristic equation 

as all distinct. Let Zi represent certain determinate functions of x of the 

form 

Zi{x)= e-^*<^'a;-"" "P^ {x) , 
where 

, , , ^^0;*^+^ , x^ , x^-'^ 

fi(^^ = ~j^i + "-'--^i^ +"'.-*+i~fc^ir + • • • • +H-i^y 

and 

.,n. 



P,W~[«.,.+ ^ + ^ + •■••]! i=l. 



when X is large, real, and positive. Then for the same values of x there exists 
a particular solution, y{x), of (1), ivhich can he developed asymptotically as 
follows: 

(a) y{x)~x'^-^AA,+ ^-\-^ + ...A,iffi{x)^^, i^l,....,n. 



and 



5g+-i + _| + J, if f Ax) ^0, where r is one 

only one of the values 1, 2, ....,«, and if at the same time m—oi.,.o^S, 
when S takes on the values — 1, 0, 1, 2. . . . 

(c) ^/(ct;)~^'"-'"-"''+i[5o+ 1^ + § + • • • .]+^""°-"'-"'[co+ ^^+....~^logx, 

if f^x) = where r is one and only one of the values 1, 2, . . . ., n, and if 
at the same time m — a^^Q—S when 8 is minus one, zero, or a positive 
integer. 

University of Michigan. 



